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The role of miRNA and miRNA biogenesis genes in the adult brain is just beginning to be 
explored. In this study we have performed a comprehensive analysis of the expression, 
genetic regulation, and co-expression of major components of the miRNA biogenesis 
pathway using human and mouse data sets and resources available on the GeneNetwork 
web site (genenetwork.org). We found a wide range of variation in expression in both 
species for key components of the pathway — Drosha, Pasha, and Dicer. Across species, 
tissues, and expression platforms all three genes are generally well-correlated. No single 
genetic locus exerts a strong and consistent influence on the expression of these key 
genes across murine brain regions. However, in mouse striatum, many members of 
the miRNA pathway are correlated — including Dicer, Droslia, Paslia, Ars2 (Srrt), Eif2cl 
(Ago!), Eif2c2 (Ago2), ZcclicTT, and Snip!. The expression of these genes may be partly 
influenced by a locus on Chromosome 9 (105.67-106.32 Mb). We explored ~1500 brain 
phenotypes available for the C57BL/6J x DBA/2J (BXD) genetic mouse population in 
order to identify miRNA biogenesis genes correlated with traits related to addiction and 
psychiatric disorders. We found a significant association between expression of Dicer 
and Drosha in several brain regions and the response to many drugs of abuse, including 
ethanol, cocaine, and methamphetamine. Expression of Dicer, Drosha, and Pasha in most 
of the brain regions explored is strongly correlated with the expression of key members 
of the dopamine system. Drosha, Pasha, and Dicer expression is also correlated with the 
expression of behavioral traits measuring depression and sensorimotor gating, impulsivity, 
and anxiety, respectively. Our study provides a global survey of the expression and 
regulation of key miRNA biogenesis genes in brain and provides preliminary support for 
the involvement of these genes and their product mlRNAs in addiction and psychiatric 
disease processes. 
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INTRODUCTION 

Small non-coding microRNAs (miRNAs) are a recently 
discovered class of post-transcriptional regulatory molecules that 
control gene expression. Genes involved in the biogenesis of these 
regulatory RNAs include Dicer (Dicerl), Drosha, and DgcrS 
(Pasha). As part of the canonical pathway, Drosha and Pasha 
process pri-miRNA into pre-miRNA, which is then transported 
into the cytoplasm via Exportin 5 (Bohnsack et al., 2004). 
Pasha identifies the complex stem-loop structure of pri-mRNA 
and recruits Drosha, resulting in cleavage and the generation 
of pre-miRNA (Denli et al, 2004). Once in the cytoplasm, 
pre-miRNA is processed into miRNA by Dicer, Argonaute (Ago) 
proteins, accessory proteins, and TRBP. The ribonuclease Dicer 
processes pre-miRNA into a ~21nt miRNA duplex as part of 
the RNA-induced silencing complex (RISC) loading complex 



and the duplex is further unwound by various Ago RNA-binding 
proteins in order to produce mature ~21nt single stranded 
miRNAs. Loading of miRNA into the RISC is facilitated by TRBP 
and target matching between mRNAs and miRNAs occurs in 
the RISC, which is composed of Ago proteins and accessory 
proteins. For an excellent review of the miRNA biogenesis path- 
way in brain, please see O'CarroU and Schaefer (2013). Several 
alternative and so-called non-canonical pathways for miRNA 
biogenesis also exist, including a Dicer-dependent pathway 
which functions independently of Drosha and Pasha (Berezikov 
et al, 2007; Babiarz et al, 2008), and a Dicer-independent 
pathway that is dependent on the actions of Drosha and Ago2 
(Cheloufi et al, 2010). Loss of function of any of these genes 
will thus have a severe impact of the miRNA landscape and the 
post-transcriptional regulation of mRNA (Babiarz et al, 201 1). 
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Many of the miRNA biogenesis genes are required very early 
in development both due to the essential role of miRNAs in 
these processes and because biogenesis genes function in multi- 
ple biological processes. For example, Dicer knock down results 
in early embryonic lethality with a concomitant lack of stem cell 
populations (Bernstein et al., 2003) likely due to the absence of 
key developmental miRNAs but also potentially due to defects 
in differentiation and maintenance of heterochromatin structure 
(Kanellopoulou et al, 2005). miRNA biogenesis genes also have 
important roles in neuronal development, survival, and main- 
tenance (Hsu et al, 2012; McLoughlin et al, 2012). Because of 
the critical nature of these genes during development, it is per- 
haps unexpected that they would also play roles in the mature 
central nervous system (CNS). However, it has been shown that 
genetic variants in miRNA processing genes have been linked to 
schizophrenia (Zhang et al., 2012; Zhou et al, 2013), depression 
(He et al, 2012), and Huntington's disease (Lee et al., 2011). In 
addition, a recent study by Tapocik et al. (2009) identified Dicer 
as one of the possible candidate genes mediating differential mor- 
phine analgesia between C57BL/6 (B6) and DBA/2J (D2) mouse 
strains. Addiction, psychiatric disorders, and neurological disease 
are influenced by genetic factors and these initial studies suggest 
that miRNA and miRNA biogenesis associated genes may play a 
role in the genetic regulation of these disorders. 

In the present study, we provide the first comprehensive anal- 
ysis of variation in miRNA biogenesis gene expression across 
both mouse and human populations and investigate the poten- 
tial behavioral and physiological consequences of this variation. 
First, we use a large family of genetically diverse strains derived 
from a cross between B6 and D2 — the BXD family — to dis- 
sect sources and potential consequences of miRNA biogenesis 
gene variation. Currently the largest and most well-characterized 
genetic family of mice, the BXD set of strains has been profiled 
for ~1500 brain and behavioral phenotypes and the whole fam- 
ily is densely genotyped. This set of ~80 strains has also been 
extensively profiled for gene expression across 12 brain regions 
and sub-regions using a variety of platforms. Finally, we com- 
pare the expression of miRNA biogenesis genes in human brain 
from normal, Alzheimer's (AD), and Huntington's (HD) disease 
cases. This wealth of molecular, genetic, and phenotypic data 
is freely available at the GeneNetwork (www.genenetwork.org) 
web site. We leverage the information available at this excellent 
functional genomics resource to investigate the genetic regula- 
tion, expression, and coexpression networks of miRNA biogenesis 
genes and to dissect their contribution to disease, addiction and 
psychiatric-related behaviors. 

METHODS 

GeneNETWORK DATABASE DESCRIPTIONS 

The BXD family includes ~80 recombinant inbred strains derived 
at three different time points from separate crosses between a 
female B6 and a male D2. The Fl progeny were subsequently 
intercrossed followed by inbreeding to fix parental genotypes at 
each locus. The first set (BXDl through 32) was created in the 
late 1970s and the second set (BXD33 through BXD42) was gen- 
erated in the early 1990s, both by Benjamin Taylor. The third set 
(BXD43 through 100) was created in the late 1990s and early 



2000s at UTHSC (Peirce et al, 2004). The following BXD brain 
expression data sets were used in this analysis: Hippocampus 
Consortium M430v2 (Jun06) RMA (Overall et al, 2009), INIA 
Amygdala Cohort Affy MoGene 1.0 ST (Marll) RMA, HBP 
Rosen Striatum M430V2 (AprOS) RMA Clean, VCU BXD PFC 
Sal M430 2.0 (Dec06) RMA, VCU BXD PFC EtOH M430 2.0 
(Dec06) RMA, VCU BXD VTA Sal M430 2.0 (Iun09) RMA, VCU 
BXD VTA EtOH M430 2.0 (Iun09) RMA, VCU BXD NA Sal 
M430 2.0 (Oct07) RMA, and VCU BXD NA EtOH M430 2.0 
(Oct07) RMA. All data sets have been standardized to a mean 
log2 expression level of 8 with a standard deviation of 2. Average 
log2 expression less than 6 is considered to be background noise. 
Detailed information for each data set is available at genenet- 
work.org. 

The following human brain data sets were also 
included,"GSE5281 Human Brain Full Liang (Jul09) RMA," 
"GSE15222 Human Brain Meyers (Apr09) Ranklnv", and 
Harvard Brain Tissue Resource Center "HBTRC-MLC Human 
Prefrontal Cortex Agilent (Junll) mlratio." For the Liang data 
set, cases include ~16 normal aged subjects and ~16 AD subjects 
(Liang et al., 2006, 2008), and gene expression was profiled for 
sis regions on the Affymetrix platform. For the Meyers data 
set, expression was profiled for temporal cortex and cortical 
tissue from 187 normal aged adults and 176 AD cases on the 
lUumina Sentrix Bead array (HumanRef-8) using lUumina's rank 
invariant transform (Webster et al., 2009). The HBTRC data 
set includes expression profiles for prefrontal cortex from ~170 
normal subjects, ~230 HD subjects, and ~400 AD subjects using 
a custom-made Agilent 44K microarray. The Harvard Brain 
dataset was contributed by Merk Pharmaceutical through the 
Sage Bionetworks Repository. Additional information for each 
data set is available at genenetwork.org. 

COVARIATION OF mlRNA BIOGENESIS GENES AND CENTRAL 
NERVOUS SYSTEM PHENOTYPES 

The BXD phenotype database was used to find the top 100 corre- 
lates of Pasha, Drosha, and Dicer expression in each brain region. 
The Spearman correlation coefficient was used to reduce the 
effect of outliers and a moderately strict criterion of significance 
of P < 0.01 was used to select the top correlates. In the case 
of multiple probe sets targeting different mRNA regions of the 
same gene, the highest expressed probe set was used. An adjusted 
p-value was calculated based on 20,000 permutations for the 
Central Nervous System (CNS) phenotype correlations to address 
the large number of phenotypes available in GeneNetwork. The 
adjusted p-value was used to sort and prioritize CNS phenotypes. 
The criterion for determining a significant correlation was an 
adjusted p-value less than 0.01 after permutation tests (Table S3) 
or a p-value less than 0.005 (Table S2). 

mlRNA BIOGENESIS GENE SELECTION AND NETWORK 
CONSTRUCTION 

Mouse genes with the following Gene Ontology (GO) 
categories — primary miRNA processing (GO:0031053), pre- 
miRNA processing (GO:0031054), and production of miRNAs 
involved in gene silencing by miRNA (GO:0035196) — were 
selected using Amigo (amigo.geneontology.org). Literature 
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searches were used to confirm, expand, and refine categories. 
Unique gene symbols were used to retrieve expression data 
fi-om each GeneNetwork data set. Probe sets targeting each 
gene were screened based on specificity (probe sets target the 
gene of interest and do not overlap polymorphisms that affect 
hybridization (see Ciobanu et al., 2010 for review), expression 
(generally greater than 8 log2 expression units), and target [probe 
sets targeting exons have less complicated expression patterns and 
network correlations compared to probe sets targeting 3' UTR 
regions, see (Mulligan et al., 2012) for review]. AH correlations, 
network graphs, and principal component analyses (PCA) were 
performed using tools and functions available in GeneNetwork. 
Briefly, PCA is implemented using custom Python code and is 
automatically computed for correlation matrices with three or 
more members. Output of the PCA function are a list of trait 
vectors representing the top principal components that can be 
used with any GeneNetwork tool, a Scree Plot, and a Factor 
Loading Plot. 

EXPRESSION OF miRNA BIOGENESIS GENES AFTER ETHANOL 
EXPOSURE IN BXD STRAINS 

We used expression data collected from the prefrontal cortex, ven- 
tral tegmental area, and nucleus accumbens of BXD strains 4 h 
following either a saline injection (the "VCU Sal" data sets) or a 
1.8 g/kg interperitoneal injection of ethanol ("VCU EtOH" data 
sets) (Kerns et al., 2005). Biological response to ethanol without 
accounting for the effect of genotype was measured by an uncor- 
rected paired f-test for each region between saline and ethanol 
treatment groups. A more sophisticated analysis (Fisher's Exact 
Test) that takes differential response to treatment due to the effect 
of genotype into account has been described previously for this 
data set (Wolen et al., 2012) and genes significantly regulated by 
ethanol treatment are listed in Table SI of that publication. We 
queried this table for alcohol responsive miRNA biogenesis genes. 

RESULTS 

EXPRESSION AND GENETIC REGULATION OF mlRNA BIOGENESIS 
GENES IN MOUSE BRAIN 

Drosha, Pasha, and Dicer expression is measured from multiple 
probe sets that target different mRNA features. Probe sets that 
target coding exons and the 3' UTR are well-expressed across 
the amygdala, prefrontal cortex, striatum, and hippocampus ( > 8) 
while probe sets targeting introns are generally not well-expressed 
(<6) and are excluded from future analysis (Table 1). An excep- 
tion occurs for the last four coding exons of PasJia ( 142298 l_at), 
which are well-expressed in the hippocampus but not in the pre- 
frontal cortex and striatum. Probe sets do not show a strong bias 
in expression measurements caused by SNPs or other variants 
overlapping the probe target sequence and thus, are appropriate 
for use in further analysis. 

Variation among the BXD family of strains is generally modest 
(<1.4) with the exception of Drosha expression in the striatum 
and hippocampus and Pasha expression in the hippocampus (~2 
fold). None of the biogenesis genes show consistent and signifi- 
cant local modulation of expression by sequence variants located 
within or near each gene locus (a cis eQTL) across multiple probe 
sets and brain regions (Figure 1). However, expression of the last 



few exons and the 3' UTR of Drosha in the hippocampus are regu- 
lated by a significant cis eQTL and expression of the last exon and 
3' UTR is regulated by a suggestive cis eQTL in the prefrontal cor- 
tex. Each biogenesis gene is genetically uncoupled from the others 
such that the expression of Pasha, Dicer, or Drosha is not modu- 
lated from the physical location of the other two genes (a trans 
eQTL) (Figure 1). These results suggest that, with the exception 
of possible 3' UTR or splice isoforms of Drosha, expression of 
miRNA biogenesis genes in the BXD population is not controlled 
by sequence variants in cis-regulatory regions but is likely the 
result of many complex genetic interactions of small effect size. 

COEXPRESSION OF PASHA, DICER, AND DROSHA 

Biogenesis genes belong to the same biological pathway and 
expression of these key genes may be coordinately regulated. 
Indeed, Pasha, Drosha, and Dicer expression in multiple brain 
regions tends to be tightly correlated (Figures Al, A2). Average 
network correlations are 0.6 (p < 0.01), 0.5 (p < 0.01), 0.45 
(p < 0.01), and 0.3 [p < 0.05) for the striatum, hippocampus, 
prefrontal cortex, and amygdala, respectively. Network complex- 
ity arises from the existence of multiple probe sets targeting dif- 
ferent mRNA regions of the same gene (Figure Al). This is typical 
for Affymetrix M430 arrays and probe sets (Mulligan et al, 2012). 
With the exception of the hippocampus, networks constructed 
from the three key biogenesis genes using the highest expressing 
probe sets tend to be positively correlated. Coexpression of these 
three genes in each brain region cannot be explained by a single 
genetic locus, yet another indication of complex network control 
of expression of miRNA biogenesis genes. 

THE STRIATAL mlRNA BIOGENESIS NETWORK 

The striatum is a key region implicated in many CNS dis- 
eases and disorders, including psychiatric disease and addic- 
tion. It is also composed largely of a single cell type — spiny 
projection neurons — making it an excellent tissue for explor- 
ing the network connectivity of miRNA biogenesis genes. Eight 
genes involved in this pathway are connected within a sin- 
gle striatal expression network by at least two nodes with 
a \r\ > 0.30 (Figurel). Pasha (DgcrS; 1455311_at), Dicerl 
(1460571_at), Drosha (1438505_s_at), Ars2 (Srrt; 1417655_at), 
and A^oi {Eif2cl; 1434331_at) are well-connected and positively 
correlated. In contrast, Ago2 {Eif2c2; 1435636_at) and Zcchcll 
(Tut4; 1437395_at) are negatively correlated with network mem- 
bers. Snipl (1459773_x_at) membership in the network is more 
complex, with both negative and positive correlations with 
other members. In general, positive and negative correlations 
are consistent with known regulatory interaction between these 
eight genes. For example, Snipl and Ars2 are accessory proteins 
involved in pri-miRNA processing that act as positive regulators 
of Drosha-mediated processing in the canonical pathway (Krol 
et al., 2010). Agol and Ago2 are accessory proteins for the RISC 
and it is less clear why Ago2 would be negatively correlated with 
network members (Krol et al, 2010). Zcchcll acts a suppressor to 
miRNA biogenesis (Heo et al, 2009). 

A tightly correlated biological network can be reduced to an 
expression signature — or principal component — using principal 
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FIGURE 1 I QTL heat map for Dicer, Drosha, and Pasha across brain 
regions from tlie BXD population. Warm colors indicate higher expression 
from the D2 allele and cooler colors indicate higher expression associated 
with inheritance of the B6 allele. Arrowheads indicate genomic location of 
each mlRNA biogenesis gene. A limited number of Drosha probe sets are 



associated with significant (hippocampus) or suggestive (prefrontal cortex) 
regulation of expression by local sequence variants — a cis eQTL. Both probe 
sets have higher expression in BXD strains that inherited the B6 allele at this 
locus. INIA_AmgCoh, Amygdala; VCUSal, Prefontal Cortex; SA, Striatum; HC, 
Hippocampus. 
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FIGURE 2 I Network correlations for miRNA biogenesis genes in the 
striatum. The network shown in the left panel was constructed using 
the highest expressing probe sets from the striatum data set IHBP 
Rosen Striatum M430V2 (Apr05) RMA Clean]. Correlation strength is 
indicated by network edge color with warm colors indicating positive 
correlations and cool colors indicating negative correlations. Common 
aliases: Eif2c1 {AgoD, Eif2c2 {Ago2}, Zcchcll {Tut4), and Ars2 {Srrt). 
Results of a principal component analysis are shown in the top right 
box. Principal components or loading factors are shown on the X-axis 



and percentage of variance explained is shown on the Y-axis. A 
suggestive eQTL on Chr9 controls the expression of principal component 
1 (lower right box). The blue line plots the LRS value by physical 
position on Chr9. The actual LRS value is shown on the left Y-axis and 
the additive effect is shown on the right Y-axis. Yellow tick marks on the 
bottom X-axis show position of SNPs and colored bars on the upper 
X-axis show positions of genes. Red and gray horizontal lines mark the 
threshold for genome-wide significant (P < 0.05) and suggestive 
thresholds based on 2000 permutations. 



component analysis (PCA). For the striatal miRNA biogene- 
sis network we used the first principal component (PCI) as a 
network signature to identify phenotypes that covary with net- 
work expression. PCI for the miRNA biogenesis network in the 
striatum explains ~50% of the variance in network coexpres- 
sion (Figure 2). PCI expression maps to a region on Chr 9 
(105.67-106.32 Mb) with a LRS of 9.6, which meets the genome- 
wide suggestive threshold (Figure 2). The top behavioral pheno- 
types correlated with the expression of PCI {p < 0.005) are all 
psychiatric- or addiction-related traits (Table SI). This includes 
anxiety (12463, 12476, 12381, 12460, 12478, 12464, 12462, 12383) 
and sensorimotor gating traits (11942) and locomotor responses 
to cocaine (10488), ethanol (11701), and morphine (11342, 
11856, 11343). Taken together these results suggest a role for 
miRNA biogenesis genes and miRNA pathways in psychiatric- 
and addiction-related phenotypes. 

COVARIATION WITH ADDICTION AND BEHAVIORAL PHENOTYPES 

Expression of each miRNA biogenesis gene is significantly cor- 
related with many central nervous system phenotypes in the 
amygdala, striatum, prefrontal cortex, and the hippocampus 
(Figures and Tables S2, S3). These brain regions are important 
for reward, decision-making, and emotional response. Alterations 
in gene expression and neuronal circuitry in these regions have 
been implicated in psychiatric disorders and addiction. To bet- 
ter assess the involvement of Drosha, Dicer, and Pasha in these 



disorders, we examined the correlation between expression of 
each biogenesis gene and the ~1500 BXD CNS phenotypes avail- 
able at the GeneNetwork web resource. As shown in Figure 3 for 
each brain region, the expression of miRNA biogenesis genes is 
highly significantly correlated with addiction-related phenotypes. 
The expression of Dicer and Drosha is correlated with cocaine, 
ethanol, and methamphetamine responses. Pasha is correlated 
with endophenotypes relevant to drug taking behaviors, such as a 
measure of impulsivity and response to a novel environment. 

The dopamine system has been strongly linked to addiction 
processes, particularly the reward aspects of addiction (Wise, 
2009) and all of the miRNA biogenesis genes exhibit strong 
correlations with various phenotypes related to the dopamine 
system across most of the brain regions surveyed. For example. 
Pasha and Dicer expression in the hippocampus is correlated with 
response to a dopamine D2-type (Drd2) and D3-type receptor 
agonist (10048) and dopamine transporter levels in the ventral 
midbrain (10282), respectively. Prefrontal cortical expression of 
Dicer and Pasha is correlated with DRD2 and dopamine trans- 
porter levels in the same tissue (10280 for Dicer; and 10267, 
10279, 10281 for Pasha). Additionally, the expression of Dicer and 
Drosha in the striatum is correlated with striatal volume (13440 
for Dicer and 11498 and 10998 for Drosha) suggesting a possible 
morphological link. 

The expression of miRNA biogenesis genes is also correlated 
with endophenotypes for psychiatric disorders. For example. 
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FIGURE 3 I CNS phenotypes strongly correlated with miRNA biogenesis 
brain gene expression. Detailed phenotype descriptions can be found in 
Table S2 and Table S3. Top correlated phenotypes were selected for 
Hippocampus, Amygdala, Prefrontal Cortex, and Striatum based on 
permutation testing (adjusted P-value < 0.01) and a network was constructed 
for each region using tools in GeneNetwork. Negative network correlations 
are indicated by cool line colors and positive correlations are indicated by 
warm line colors. Correlation threshold is r = 0.5. At the center of each 
network, Pahsa (DgcrSj, Drosha, and Dicer nodes are shaded in blue. Yellow, 
purple, green, and red nodes indicate addiction-, psychiatric-, dopaminergic-. 



and morphology-related traits, respectively. In each brain region, with the 
exception of the amygdala, many behavioral phenotypes related to addiction 
and psychiatric disease are highly correlated with the expression of miRNA 
biogenesis genes. Importantly, dopaminergic phenotypes and morphological 
phenotypes in key brain regions that play a role in addiction and psychiatric 
disease processes are also correlated with behavioral phenotypes and 
expression of miRNA biogenesis genes. This indicates an overlap between 
miRNA biogenesis gene expression and developmental processes and 
biological pathways involved in addiction and psychiatric disease-related 
phenotypes. 



Drosha expression in the hippocampus is correlated with mea- 
sures of depression (12554). Prefrontal cortical expression of 
Pasha is correlated with a measure of impulsivity (12731) 
and Dicer levels in the striatum are correlated with anxiety 



measurements (12477, 12478, 12451, 12476). Additionally, stri- 
atal levels of Drosha are correlated with pre-pulse inhibition of 
the acoustic startle response (11942, 11684, 11941) which is a 
measure of sensorimotor gating — a process disrupted in many 
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schizophrenic patients. Our results show that the expression of 
miRNA biogenesis genes is directly related to numerous pheno- 
types associated with addiction and psychiatric disease as well as 
other molecular pathways linked to these phenotypes. 

EXPRESSION OF miRNA BIOGENESIS GENES AFTER ETHANOL 
EXPOSURE IN BXD STRAINS 

To directly assess whether exposure to ethanol alters the expres- 
sion or genetic regulation of miRNA biogenesis genes, we queried 
expression data sets collected from the prefrontal cortex, ven- 
tral tegmental area, and nucleus accumbens of BXD strains four 
hours following either a saline injection or a 1.8 g/kg inter- 
peritoneal injection of ethanol (Kerns et al., 2005). A signifi- 
cant expression difference between treatment groups — without 
accounting for genotype — was detected by paired f-test (uncor- 
rected, p < 0.05) in the nucleus accumbens for Pasha probe 
sets targeting the 3' UTR (1435439_at) and an intron with low 
expression (1446169_at). A recent analysis allowing for differen- 
tial response to ethanol due to genetic variation among the BXD 
strains identified Ago2 (1426366_at) and Pasha ( 14553 ll_at) 
probe sets in the prefrontal cortex and Ars2 (1417655_a_at) 
and Ago2 (1426366_at) probe sets in the nucleus accumbens as 
being regulated by ethanol (Wolen et al., 2012). Significant cis 
modulation of expression (LRS = 16.2) for a Pasha probe set 
( 14553 ll_at) targeting the distal 3' UTR was also detected after 
ethanol treatment — but was absent in saline injected controls — 
in the nucleus accumbens. Exposure to alcohol significantly 
alters the genetic modulation oi Dicer (1427941_at; Chr9@47.04 
Mb; LRS = 16.97) in the ventral tegmental area and both 
Dicer (1460571_at; Chr9@69.81 Mb; LRS = 15.6) and Drosha 
(1438505_s_at, Chr9@65.31 Mb; LRS = 15.7) in the prefrontal 
cortex. As observed in other naive BXD brain expression data sets, 
no significant modulation of expression was observed in control 
BXD strains. These results suggest that miRNA biogenesis gene 
expression can be modulated by exposure to alcohol. 

EXPRESSION OF mlRNA BIOGENESIS GENES IN HUMAN BRAIN 

Expression levels of miRNA biogenesis genes — DICERl, DGCR8 
(PASHA), and DROSHA (RNASEN or RNASE3L) show a remark- 
able level of variation across normal aged human brain (~ 1.3-68 
fold difference) using different platforms, populations, and probe 
sets (Table S4). Human populations are incredibly polymorphic 
and some of this variation reflects true expression differences as 
well as technical artifacts due to genetic variants that overlap the 
probe target sequence and alter binding. In many of the human 
data sets, there is considerable covariation of PASHA, DROSHA, 
und DICERl (Figure A3). 

To assess the involvement of miRNAs and miRNA biogenesis 
genes in human neurological diseases, we explored the expres- 
sion of these genes across several brain regions from control 
and Alzheimer's (AD) and Huntington's (HD) disease cases. AD 
is the major cause of adult dementia and involves early loss of 
neurons — especially cholinergic — in the cortex and hippocam- 
pus. AD is also associated with a wide range of psychiatric 
symptoms (Koppel et al., 2012). HD involves the selective loss 
of striatal populations of medium spiny projection neurons and 
concomitant alterations in mood — often depression — and motor 



coordination. Intriguingly, we found highly statistically signifi- 
cant expression differences for miRNA biogenesis genes between 
control and AD or HD disease cases (Figure 4 and Table S4) based 
on uncorrected f-tests between control and disease cases. For HD 
cases compared to controls, most PASHA {DGCR8) probe sets 
and DICERl have higher expression in the prefrontal cortex while 
DROSHA expression is decreased. For AD cases compared to 
controls, there is generally significantly (p < 0.05) higher expres- 
sion of PASHA in prefrontal cortex, neocortex, and entorhinal 
cortex and significantly lower expression in the hippocampus, 
medial temporal gyrus, cingulate cortex, and superior frontal 
gyrus. Significantly lower expression of most DICER probe sets 
(212888_at, 213229_at) is detected in AD cases compared to con- 
trols and significantly lower expression of DROSHA in AD cases 
compared to controls is also observed in most data sets (HBTR 
and Meyers) (Figure 4). Taken together, this survey of miRNA 
biogenesis gene expression across brains regions from normal 
or disease cases suggests an association between expression of 
these genes and AD and HD. However, the mechanism underly- 
ing this relationship and whether or not the observed alterations 
in miRNA biogenesis gene expression are causal or downstream 
of the disease state are not known. 

DISCUSSION 

In this study, we used bioinformatic and genetic analyses of 
miRNA biogenesis genes to show that microRNAs could be crit- 
ical modulators of pathways that underlie differential responses 
to ethanol or drugs of abuse and differential indicators of psychi- 
atric disease. Our initial analysis examined the expression of these 
genes and showed that there were both strain-specific and brain 
region-specific differences. This wide range of variation in expres- 
sion was found in both mouse and human populations. Despite 
this variation, the key components of the pathway — Drosha, 
Pasha, and Dicer — are generally strongly correlated across brain 
regions. Subsequent analyses were used to assess whether these 
expression differences could impact phenotypes. For example. 
Dicer, Drosha, and Pasha are correlated with key CNS pheno- 
types related to addiction and psychiatric disorders. Additionally, 
exposure to alcohol alters the expression and genetic regulation of 
these genes suggesting that differential responses between strains 
as well as a priori expression differences across strains could be 
critical in mediating phenotypic differences. 

We observed a large level of variation in the expression of 
Drosha and Pasha in the striatum and hippocampus. Using the 
BXD genetic mouse resource we can rule out technical sources of 
variation and focus on possible mechanisms of genetic regulation. 
Interestingly, none of these genes are consistently modulated by 
local sequence variants that influence expression, at least in brain. 
It also appears that there is no single trans modulatory locus 
that controls a major component of their expression. Even in 
the striatum, where there is a relatively homogenous population 
of cells and strong correlations among eight miRNA biogene- 
sis pathway members, the best single locus modulating network 
expression (Chr 9 from 105.67 to 106.32 Mb) explains less than 
50% of the joint variation in expression. Candidate genes within 
this region — Atp2cl, Tlr9, and Alasl — may have a broad, yet sub- 
tle, influence on miRNA biogenesis pathways. Indeed, many genes 
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FIGURE 4 I Expression of miRNA biogenesis genes across human 
cases. Expression of multiple DICER, DROSHA, and PASHA probe sets 
from different regions of normal and diseased human brain tissue. (A) 
Expression levels are shown for the three mlRNA biogenesis genes in 
the prefrontal cortex of normal or Alzheimer (AD) or Huntington Disease 
(HD) cases (HBTRC data). The average expression ratio of each gene in 
normal, AD, and HD cases compared to a common reference sample is 
shown on the Y-axis. (B) Average log2 expression (Y-axis) of DICER. 
PASHA, and DROSHA is shown for the Neocortex of normal or AD 



and steps are involved in the production of miRNA and the 
process is tightly regulated. A major mutation in this pathway 
severely compromises organism function very early in develop- 
ment. The BXD genetic reference population is not segregating 
for a major mutation in key miRNA biogenesis genes, but proves 
useful to assess how subtle variation influences their expression 
and behavioral phenotypes. 

Across brain regions and platforms, expression of Pasha, Dicer, 
and Drosha significantly covaries with the expression of mul- 
tiple addiction phenotypes and psychiatric-related traits. Both 
Dicer and Drosha were correlated with the response to cocaine, 
ethanol, and methamphetamine. This suggests that these genes 
are related to the behavioral response to a wide-range of sub- 
stances and are more related to addiction, per se, rather than 
to any specific substances. In addition, the category or type of 
response shows considerable consistency. For example, locomo- 
tor activities are correlated with expression of all three enzymes 



patients. Expression data was standardized to have an average log2 
expression of 8 with a standard deviation of ± 2. (C) Average log2 
expression of all three genes, as assayed by different probe sets, is 
shown on the Y-axis for 6 different brain regions from normal and AD 
cases. Expression data was standardized to have an average log2 
expression of 8 with a standard deviation of ± 2. Significant expression 
differences are often observed for at least one probe set between 
normal and Alzheimer or Huntington Disease cases. Significance was 
determined by Student's two-tailed f-test. 



in multiple brain regions and to multiple substances. This sug- 
gests that certain responses may have a stronger relationship 
to the expression of the miRNA biogenesis genes. In addition, 
the expression and genetic modulation of Pasha in the nucleus 
accumbens is significantly altered after ethanol exposure in BXD 
strains. The same exposure to alcohol results in trans modula- 
tion of Dicer and Drosha by loci on Chr9 at 47.04, 65.31, and 
69.81 Mb in the ventral tegmental area and the prefrontal cortex. 
A large region on Chr9 (~37 Mb to 99 Mb) has previously been 
associated with alcohol preference (Alpq3, Ap5q, Alpql) and con- 
sumption {Etohc3) in B6D2-derived populations, including BXD 
strains (Phillips et al, 1994, 1998; Belknap et al., 1997; Tarantino 
et al, 1998; Belknap and Atkins, 2001; Kerns et al, 2005; Bice 
et al, 2006, 2008, 2009; Mulligan et al, 2006; Weng et al., 2009). 
The alcohol QTL region on Chr 9 in mouse is syntenic to a region 
on human Chrl5 (Ehlers et al, 2010) and association studies 
in human populations have identified alcohol dependence QTLs 
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within this region of synteny at 50-53 Mb (Ehlers et al., 2004). 
These results indicate that alcohol exposure changes the expres- 
sion and genetic regulation of miRNA biogenesis genes through 
an unknown mechanism. This mechanism may influence human 
alcohol dependence as well. 

Intriguingly, all three key biogenesis genes are also correlated 
with behavioral measures of depression, anxiety, and schizophre- 
nia. There is high comorbidity between psychiatric disease and 
addiction (Swendsen et al., 2010) suggesting that these disorders 
involve similar brain pathways. Although the exact mechanism 
explaining the correlation cannot be addressed in this analysis, 
the strong covariation supports the involvement of miRNA bio- 
genesis genes in the expression of these phenotypes and warrants 
deeper scrutiny. 

The striatum is a brain region important in motivation, deci- 
sion, and reward as well as addiction, psychiatric, and mood 
disorders. We observed a strong correlation between expression 
of key members of the dopamine pathway and miRNA biogenesis 
gene expression in most brain regions. The strongest expression 
signature from eight correlated biogenesis genes in the stria- 
tum was strongly and significantly associated with locomotor 
responses to cocaine, ethanol, and morphine as well as anxiety 
and sensorimotor gating traits. miRNA and miRNA biogenesis 
genes have been previously linked to maintenance of the reward 
pathway and psychiatric disease and may have an exquisite inter- 
action with neuronal function in the striatum. For example, loss 
of Dicer in dopaminergic neurons and Drd2 expressing medium 
spiny neurons leads to death of those neuronal populations (Kim 
et al., 2007; Davis et al, 2008). Loss of A^o2 in adult medium spiny 
neurons expressing Drd2 rendered mice insensitive to cocaine and 
they failed to develop dependence (Schaefer et al, 2010). Alcohol 
consumption in adult mice leads to an approximately 2-fold 
increase in the level of miR-9 in the striatum followed by a reduc- 
tion in expression of a subset of mRNA isoforms that contain 
miR-9 targets in the 3' UTR and code for a large conductance cal- 
cium and voltage gated potassium channel (Pietrzykowski et al, 
2008). In human populations, microdeletions that include the 
DGCR8 gene (PASHA) cause DiGeorge syndrome — a multispec- 
trum disorder associated with a high risk of schizophrenia. We are 
still exploring the role of miRNA in adult brain but it may be that 
this class of regulatory molecules and their product miRNAs are 
functionally more critical in certain regions and cell types then 
others. 

Our analysis of miRNA biogenesis gene expression and reg- 
ulation was primarily based on extensive data collected from 
the BXD genetic mouse population. The majority of our find- 
ings in mouse also extend to human brain tissue where we 
observed an even greater level of expression variation and 
coexpression of PASHA, DICER, and DROSHA. These genes are 



highly polymorphic in humans and without genotyping we can- 
not eliminate the possibility that expression variation is caused 
by technical artifacts related to probe design and probe target 
sequence. Interestingly, we also found highly statistically signifi- 
cant expression differences for miRNA biogenesis genes between 
control and Alzheimer's (AD) and Huntington's (HD) disease 
cases (Figure 4). It is highly unlikely that SNPs and other variants 
overlapping probe sets could cause this pattern of miRNA bio- 
genesis gene expression between control and disease cases. These 
differences also suggest that miRNA processes play some role in 
the progression of HD and AD, either as drivers or as downstream 
consequences. Indeed, down-regulation of several miRNAs has 
been documented in HD (Johnson et al., 2008; Packer et al., 
2008), which features early loss of medium spiny neuronal cell 
populations in the striatum. However, little is known about the 
involvement of miRNAs and miRNA biogenesis genes in AD. 
Our results provide additional support for the relevance of these 
genes in adult human brain and are further impetus for examin- 
ing their role in neurological disease, addiction, and psychiatric 
disorders. 

Taken together, our study provides evidence for the role 
of miRNA biogenesis genes and pathways in addiction- and 
psychiatric-related disease processes and suggests that miRNA 
pathways may be perturbed in human brain degenerative disor- 
ders such as HD and AD. Finally, our analysis provides additional 
support for an exquisite and critical role between miRNA and 
miRNA biogenesis and the function of the striatum. Most of 
our analysis is based upon correlations, which do not neces- 
sarily indicate direct causal relationships. The next step will be 
to establish causal links between the expression of miRNA and 
key biogenesis genes and these important CNS phenotypes and 
diseases. 
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FIGURE Al I Network correlations for all Dicer, Pasha, and Drosha {DgcrS) probe sets. Correlations shown are for all probe sets with average log 2 
expression greater than 7 from each brain region data set. 
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FIGURE A2 I Network correlations for Dicer, Pasha, and 
Drosha {DgcrS). Correlations shown are for the highest 
expressing probe sets from each brain region data set. The 



exception is the hippocampus for which the highest pairwise 
correlations among all probe sets were selected (1422981_at, 
1460571_at, 1438440_at). 
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FIGURE A3 1 Coexpression of key miRNA biogenesis genes in normal 
hiuman brain data sets. Network correlations are especially strong across 
multiple brain expression data sets from the study by Liang and colleagues 
(2006, 2008). Network connections are less robust in the data sets 


contributed by Meyers and colleagues (Webster et al., 2009). Because of the 
large number of variants and differences in probe design across platforms, 
probe sets were chosen based on expression and overall network 
connectivity. 
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